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Abstract. The stellar- wind structure in high-mass X-ray 
binaries (HMXBs) is investigated through modelling of 
their ultraviolet (UV) resonance lines. For the OB super- 
giants in two systems, Vela X-1 and 4U1700— 37, high- 
resolution UV spectra are available; for Cyg X-1, SMC 
X-1, and LMC X-4 low-resolution spectra are used. In or- 
der to account for the non-monotonic velocity structure of 
the stellar wind, a modified version of the Sobolev Exact 
Integration (SEI) method by Lamers et al. (1987) is ap- 
phed. 

The orbital modulation of the UV resonance lines provides 

information on the size of the Stromgren zone surrounding 

^^ I the X-ray source. The amplitude of the observed orbital 

PsJ ■ modulation (known as the Hatchett-McCray effect), how- 

,— H [ ever, also depends on the density- and velocity structure 

, ^ • of the ambient wind. Model profiles are presented that 

illustrate the effect on the appearance of the HM effect 

by varying stellar- wind parameters. The q parameter of 

, I Hatchett & McCray (1977), as well as other parameters 

^ ■ describing the supergiant's wind structure, are derived for 

[■^^ the 5 systems. The X-ray luminosity needed to create the 

^^D ■ observed size of the Stromgren zone is consistent with 

the observed X-ray flux. The derived wind parameters are 

compared to those determined in single OB supergiants of 

similar spectral type. 

f^ . Our models naturally explain the observed absence of the 

■^ ; HM effect in 4U1700-37. The orbital modulation in Vela 

O ■ X-1 indicates that besides the Stromgren zone other struc- 

' ' tures are present in the stellar wind (such as a photo- 

»~ ionization wake). The ratio of the wind velocity and the 

"^ ■ escape velocity is found to be lower in OB supergiants in 

Ci • HMXBs than in single OB supergiants of the same effec- 

^ ' tive temperature. 

• I— I 

^\ Key words: Line: profiles — binaries: close — circum- 
t^ ■ stellar matter — Stars: early-type — Stars: mass-loss — 
■ " " ' Ultraviolet: stars 
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1. Introduction 

In a high-mass X-ray binary (HMXB) a massive, early- 
type star transfers mass onto a compact companion, a 
neutron star or a black hole. The potential energy asso- 
ciated with the accreted matter is efficiently converted 
into X-rays (Davidson & Ostriker 1973). The mass trans- 
fer can take place in two different ways: (i) via the massive 
star's dense stellar wind which is partly intercepted by the 
strong gravitational field of the compact companion; or 
(ii) via a flow of matter towards the compact star through 
the inner Lagrangian point (Roche- lobe overflow). In the 
latter case an accretion disk is expected to be present in 
the system, as well as a rapidly spinning X-ray source. 
The massive star is either an OB supergiant with a dense 
radiation-driven wind or a Be-type star characterised by 
strong Ha emission arising from a slowly outflowing, dense 
equatorial disk. For reviews on this subject we refer to 
Lewin et al. (1995), van Paradijs (1998) and Kaper (2001). 
The accretion-induced X-ray luminosity gives a direct 
measure of the wind density and velocity at the orbit 
of the X-ray source. In this sense, the compact object 
acts as a probe in the stellar wind. In OB-supergiant sys- 
tems, the X-ray source strongly affects the supergiant's 
radiation-driven wind. The X-ray source ionizes the sur- 
rounding wind regions creating an extended (Stromgren) 
zone of strong ionization that trails the X-ray source in 
its orbit. This causes the orbital modulation of ultra- 
violet (UV) resonance lines (Hatchett & McCray 1977; 
Kaper et al. 1993). Hatchett & McCray predicted this ef- 
fect to be observable in UV resonance lines of the HMXB 
HD153919/4U1700-37. The HM-effect has, however, not 
been detected in spectra obtained with the International 
Ultraviolet Explorer (lUE), as shown by Dupree et al. 
(1978). Although Hammerschlag-Hensberge et al. (1990) 
reported variations with orbital phase in some subordinate 
lines, Kaper et al. (1990) reported that these variations are 
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caused by variable Raman-scattered emission lines and are 
not due to the HM-effect (see also Kaper et al. 1993). 

Dupree et al. (1980) did detect the HM-effect in 
lUE spectra of HD77581/Vela X-1. The B-supergiant's 
wind profiles are less saturated while the Stromgren 
zone is expected to be larger than in the case of 
HD153919/4U1700-37. A detailed study of the orbital 
modulation of the UV resonance lines of HD77581/Vela 
X-1 indicated that the velocity and density structure of 
the stellar wind cannot be a monotonically rising function 
with distance from the star (Kaper et al. 1993), a con- 
clusion that has been derived from observations of single 
early- type stars as well (Lucy 1982). Additional absorp- 
tion (e.g. due to material trailing the X-ray source in its 
orbit) appears in the line profiles at late orbital phases 
(Sadakane et al. 1985; Kaper et al. 1994). 

HD153919 and HD77581 are the only OB-supergiants 
in HMXBs bright enough in the UV to have been ob- 
served with sufficient signal-to-noise in the high-resolution 
mode of lUE. Kallman et al. (1987) and Payne & Coe 
(1987) studied the HD77581/Vela X-1 system in low- 
resolution mode and trailing mode, respectively, to search 
for the signature of rapid modulation of the Stromgren 
zone as a reaction to the periodically varying X-ray emis- 
sion from the pulsar. Although these early experiments 
with lUE were unsuccessful, the periodic modulation was 
recently detected using the Faint Object Spectrograph on- 
board the Hubble Space Telescope (HST) by Boroson et al. 
(1996). Treves et al. (1980) detected the orbital modula- 
tion of UV resonance lines in low-resolution lUE spectra 
of HD226868/Cyg X-1. This source is too faint to ob- 
tain useful high-resolution lUE spectra (Davis & Hart- 
mann 1983). Bonnet-Bidaud et al. (1981) and van der 
Klis et al. (1982) discussed the appearance of the HM- 
effect in low-resolution lUE spectra of Sk 160/SMC X-1 
and Sk-Ph/LMC X-4. Two high-resolution UV spectra of 
Sk 160/SMC X-1 confirmed the orbital modulation of the 
Si IV resonance doublet (Hammerschlag-Hensberge et al. 
1984). Vrtilek et al. (1997) and Boroson et al. (1999) dis- 
cussed observations of LMC X-4, including high-resolution 
HST/GHRS spectra, in terms of a shadow wind. In this 
system the X-ray flux is so strong that only in the X-ray 
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Table 1. System parameters of the HMXBs in our sample. References are: a) Walborn (1973); b) Chevalier & Ilovaisky 
(1977); c) Webster et al. (1972); d) Conti (1978); e) Bregman et al. (1973); f) Osmer (1973); g) Sadakane et al. (1985); 
h) Heap & Corcoran (1992); i) Bolton (1975); j) Levine et al. (1991); k) Rappaport & Joss (1983); 1) Nagase (1989); 
m) Kemp (1977); n) van Paradijs & Kuiper (1984); o) Deeter et al. (1987); p) Haberl et al. (1989); q) Bolton & Hcrbst 
(1976); r) Kaper (1998). 
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shadow behind the supergiant companion a normal stellar 
wind can develop (Blondin 1994). This represents the most 
extreme case of the HM effect; in practice, only Roche- 
lobe overflow systems will include a shadow wind. Recent 
HST/STIS observations by Kaper et al. (in preparation) 
confirm the presence of a shadow wind in LMC X-4, as 
well as a photo-ionization wake at the (leading) interface 
between the shadow wind and the X-ray ionization zone. 
We have not included the peculiar 07 Ill/black-hole can- 
didate system #32/LMC X-1 (Hutchings et al. 1983, 1987; 
Cowley et al. 1995) because the system is embedded in a 
nebula with strong (UV) emission lines. 

Although detailed observational studies of the HM 
effect in individual systems exist, a quantitative analy- 
sis based on state-of-the-art stellar-wind models is lack- 
ing. In the past two decades significant improvement has 
been made in modelling the stellar-wind profiles of (sin- 
gle) OB-type stars (Groenewegen et al. 1989; Groenewe- 
gen & Lamers 1989, 1991; Haser et al. 1998). We present 
a quantitative analysis of the UV spectral variability in 
the five HMXBs with OB-supergiant companion observed 
with lUE, comparing the complete set of lUE spectra. 
The model spectra are obtained using a modified version 
of the Sobolev Exact Integration (SEI) method introduced 
by Lamers et al. (1987). Our new method allows to take 
into account the non-monotonic wind structure observed 
in many (single) OB-type stars by adding "turbulence", 
as well as an extended Stromgren zone around the X-ray 
source. With this analysis, fundamental parameters of the 
system are derived. 

An overview of the data and spectral line variability 
is given in Sect. 2. X-ray eclipse spectra, UV continuum 
lightcurves and additional photometry are described in 
Appendices A-C, and the applied variability and error 
analysis method is described in Appendix D. Sect. 3 in- 
troduces the radiation transfer code "SEI" (Lamers et al. 
1987) and the modifications that we implemented to be 
able to compute line profiles for HMXBs. Sect. 4 describes 
our attempts to model the UV line variability observed in 
HD77581/Vela X-1 and HD153919/4U1700-37. The de- 



rived terminal velocities, ionization fractions and sizes of 
the ionization zones are discussed in Sect. 5. 

2. lUE observations of High-Mass X-ray Binaries 

Our sample comprises five HMXBs; only for 
HD77581/Vela X-1 and HD153919/4U1700-37 high- 
resolution lUE spectra will be discussed. The system 
parameters (Table 1) are collected from the literature. 
Vela X-1 is an X-ray pulsar, identifying the compact 
object as a neutron star. For 4U1700— 37 no X-ray 
pulsations have been detected, but the X-ray spectrum 
suggests that the X-ray source is a neutron star (White 
& Marshall 1983; Reynolds et al. 1999). Both HMXBs 
are wind-fed systems (Kaper 1998) and have relatively 
modest X-ray luminosities. Cyg X-1 is a well-known 
black-hole candidate. The other two sources, SMC X-1 
and LMC X-4 are located in the Magellanic Clouds. 
The low metallicity in the Magellanic Clouds and the 
occcurence of Roche- lobe overflow due to their tight orbits 
account for their high X-ray luminosities. Sk-Ph/LMC 
X-4 contains a (sub)giant rather than a supergiant. 
The short pulse period of both SMC X-1 and LMC 
X-4 suggests that these neutron stars are surrounded 
by an accretion disk. Periods of 60 and 30 days in the 
lightcurves of respectively SMC X-1 (Wojdowski et al. 
1998) and LMC X-4 (Heemskerk & van Paradijs 1989) 
are interpreted as the precession period of this disk. 

The observation logs of the lUE spectra are listed in 
Table 2. Some spectra were included in previous studies 
(Dupree et al. 1978, 1980; Treves et al. 1980; Bonnet- 
Bidaud et al. 1981; van der Klis et al. 1982; Sadakane 
et al. 1985; Hammerschlag-Hensberge et al. 1990; Kaper 
et al. 1993), but many are used here for the first time. 
The low-resolution spectra (HDE226868/Cyg X-1, Sk 
160/SMC X-1 and Sk-Ph/LMC X-4) are sampled on a 
grid of 1 A per point. The lUEDR software package (Gid- 
dings 1983) provided by STARLINK was used for data 
reduction. We refer to Kaper et al. (1993) for details 
on the data reduction of the high-resolution spectra of 
HD77581/Vela X-1 and HD153919/4U1700-37, sampled 
on a grid of 0.1 A per point. The determination and 
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Table 2. SWP number, 
spectra. Spectra used to 
parameters: 
HDE226868/Cyg X-1 
Sk 160/SMC X-1 
Sk-Ph/LMC X-4 
HD77581/Vela X-1 
HD153919/4U1700-37 



Reduced Julian Day (RJD) = JD— 2 440000 and mid-exposure orbital phase (j) of the lUE 
construct </) ~ and ~ 0.5 averages are indicated by a suffix o and -I-, respectively. Orbital 
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subtraction of the inter-order background light has been 
improved for the high-resolution lUE final archive spec- 



tra in the INES system (http://incs.vilspa.esa.es) using 
NEWSIES (Gonzalez-Riestra et al. 2000), which may be 
important for the N V resonance line around 1240 A. 

The exposure times are typically 45, 37, 45, 
150 and 30 minutes for HDE226868/Cyg X-1, Sk 
160/SMC X-1, Sk-Ph/LMC X-4, HD77581/Vela X-1 and 
HD153919/4U1700— 37, respectively, corresponding to a 
spread in orbital phase of 0.006, 0.007, 0.022, 0.012 and 
0.006 per spectrum, respectively. Average X-ray eclipse 
{4> = 0) spectra, which should be representative of the 
"undisturbed" stellar wind of the OB companion, are pre- 
sented in Appendix A. The continuum variability (which 
becomes apparent when normalising the spectra), known 



as ellipsoidal variations, is due to the "pear-like" shape of 
the supergiant filling its Roche lobe. These variations are 
discussed in Appendices B & C. The method we used for 
our error and variability analysis is based on the calcula- 
tion of covariances and is described in Appendix D. 

The strongest variability is detected in the ultravi- 
olet ^S-^P" resonance doublets that are formed in the 
stellar wind; e.g. C iv (Aq = 1548.20 & 1550.774 A, 
Aw = 498 km s'^). Si iv (Aq = 1393.755 & 1402.77 A, 
Av = 1927 km s^^), and N v (Ao = 1238.821 & 1242.804 
A, Av — 961 km s^^). The orbital modulation of the mean 
flux in these lines is shown in Figs. 1 (HDE226868/Cyg 
X-1, Sk 160/SMC X-1 and Sk-Ph/LMC X-4) and 2 
(HD77581/Vela X-1 and HD153919/4U1700-37). The er- 
ror bars follow from the covariance-based procedure. In 
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the luminous (Rochc-Lobc overflow) X-ray sourees Cyg X- 
1, SMC X-1 and LMC X-4, the X-rays ionize most of the 
stellar wind, leaving only a shadow wind unaffected. In the 
fainter (wind-fed) X-ray sources Vela X-1 and 4U1700— 37 
the stellar wind is much less disturbed. 



2.1. Cyg X-1 

In HDE226868/Cyg X-1 (Fig. 1) all three lines show clear 
modulation of the strength of the P-Cygni absorption with 
orbital phase due to the HM-effect. For C iv and Si iv this 
was already observed by Treves et al. (1980). The N v 
line shows the HM-effect too, although less pronounced. 
The modulation is very strong and not confined to orbital 
phases near — 0.5: apparently, a large fraction of the 
stellar wind is strongly ionized by the X-ray source, sug- 
gesting the presence of a shadow wind in this system. The 
low resolution does not allow a detailed study of the wind 
velocity structure, but it is clear that the modulation of 
the absorption is visible up to the highest wind velocities. 

2.2. LMC X-4 

In Sk-Ph/LMC X-4 (Fig. 1) the HM-effect is seen in aU 
three lines, in agreement with van der Klis et al. (1982). 
The N V profile is completely absent at (/> = 0.5, whereas 
part of the C iv line remains visible. High-resolution HST 
GHRS/STIS spectra (Boroson et al. 2001) show that the 
remaining C iv absorption at = 0.5 is mostly due to 
the intervening interstellar medium. These HST data also 
show some N v absorption to persist around = 0.5, 
which is probably photospheric. The strength of the P- 
Cygni lines as a function of orbital phase indicate that the 
stellar wind is confined to the X-ray shadow behind the 
OB star. The orbital modulation of the C iv and N v pro- 
files demonstrates a complexity beyond that of a smooth 
single wave, indicating a more complex mass flow. This 
has been confirmed with recent HST GHRS/STIS spectra 
(Boroson ct al. 1999; Kaper et al. in preparation). 

2.3. SMC X-1 

In Sk 160/SMC X-1 (Fig. 1) the HM-effect is seen in aU 
three lines (see also van der Klis et al. 1982). The shape 
and amplitude of the orbital modulation of the mean flux 
suggests that absorption is present in a tight orbital phase 
range around = only, i.e. a shadow wind. This is 
consistent with the high X-ray luminosity of SMC X-1 
(Hutchings 1974; van der Khs et al. 1982). 

2.1 Vela X-1 

In HD77581/Vela X-1 (Fig. 2) the HM-effect is observed 
in the C IV and Si iv fines (Dupree ct al. 1980) and the N v 
line (see Kaper et al. 1993). The severely saturated C iv 
line shows less orbital modulation than the just saturated 



Si IV line; remarkably, the N v line shows the inverse be- 
haviour. The variability in the UV lines of HD77581/Vela 
X-1 will be studied in more detail in Sect. 4.1. 



2.5. 4U 1700- 37 

In HD153919/4U1700-37 (Fig. 2) the P-Cygni fines are 
very strong and saturated. The N v line does not show 
any significant variability with orbital phase. The C iv 
and Si iv absorption lines are possibly slightly weaker very 
near (p = 0.5. This would indicate a very small Stromgren 
zone, if the variability is due to the HM-effect. The vari- 
ability in the C iv line results from changes in the blue 
edge of the P-Cygni profile, while the variability in the 
Si IV line originates at wind velocities of about —900 km 
s~^. However, the UV spectrum is affected by variable Ra- 
man scattered emission lines which show a similar orbital- 
phase dependence (Kaper et al. 1990, 1993). 

3. Modelling the UV resonance lines in HMXBs 

3.1. Adapting the SEI code 

The high-mass X-ray binary is described as a spheroidal 
early-type star with a spherically symmetric radially out- 
fiowing wind, and an orbiting X-ray point source. The 
used coordinate system is based on the impact parameter 
p and the line-of-sight parameter z, in units of the stellar 
radius and with the center of the early- type star at (0, 0): 
X = {p,z) & a; = \x\ , (1) 

and with the spectrograph at (0, — oo). The velocity v{x) of 
the radiation driven stellar wind is normalised to w(oo) = 
1 and parameterised as (Lamers ct al. 1987): 

v{x) = vo + {1 ~ vo){l - -r (2) 

X 

where vq is the velocity at the base of the wind (a; = 1 
and 7 > 0): 

Wo = w(l) ~ O.Olw(oo) , (3) 

which is of the order of the sound speed in the stellar wind. 
The X-ray point source at position xx affects the nor- 
mal abundance q'{x) of a certain ion in the wind. The 
degree at which the wind is disturbed is set by the param- 
eter q{x) (Hatchett & McCray 1977): 



q{x) 



v{x) 



v{xx) {x~ xxY 



(4) 



Approximately, at all points on a surface of constant q, 
the X-rays remove the same fraction r]{q) of the ion under 
consideration (Tarter et al. 1969). This leaves only a 1 — 
ri{q) contribution of this ion to the source function S{x) 
and optical depth t{x). In our computations we adopt a 
region enclosed by a surface set by the value of q, with 

"CM) - { 2 SL"'' " ' (5) 

A negative value for rj enhances the abundance of the ion 
within the Stromgren zone. The X-ray point source is lo- 
cated inside a closed surface of constant q if 

q < (^critical (6) 
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Fig. 1. Resonance lines of C iv (left), Si iv (middle) and N v (right) for HDE226868/Cyg X-1 (top), Sk 160/SMC X-1 
(middle) and Sk-Ph/LMC X-4 (bottom). Each upper panel displays the average spectrum near = (drawn) and 
(/) = 0.5 (dotted). The short-dashed vertical line marks the restwavelength, while the two long-dashed lines mark the 
boundaries of the region in which the mean flux is determined and plotted in the lower panel against orbital phase. 
The dotted horizontal line in each lower panel represents the mean flux level in the = average spectrum. 
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Fig. 2. Same as Fig. 1, but now for HD77581/Vela X-1 (top) and HD153919/4U1700-37 (bottom). 



with 

^critical = lim q{x) = — . (7) 

X -> oo ""(^^x) 

p = 

Because of the singularity in xx when calculating q{x), 
we excluded a small spherical region around the X-ray 
point source from the grid with a radius of one millionth 
of a stellar radius, which is of the order of the radius of 
a neutron star. As an example surfaces of constant q are 
shown (Fig. 3), for a velocity law with 

7 = 1 & vo^ 0.01 (8) 

and putting the X-ray point source at 

XX = 1.5, (9) 

which corresponds to 



The escape probability method (Castor 1970) makes 
a local approximation for the source function S^{x) at 
comoving frequency v, namely the Sobolev approximation 



S^{x) 



I3,{x)lt + eB,{x) 



(11) 



Here, (3 is the escape probability of a line photon; j3c the 
penetration probability of a continuum photon; e the ratio 
of coUisional over radiative de-excitations; B^ the Planck 
function at frequency v] and /* the stellar photospheric 
intensity at frequency v. Due to the relatively low den- 
sity of in the stellar wind, the contribution of collisional 
(de)excitations to the source function of the UV resonance 
lines can be neglected: 



V(xx) = 0.34 & (Zcritical = 2.941 . 



(10) e = 



(12) 
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Fig. 3. Sketch of the used coordinate system {p, z) indi- 
cating the location of the early-type star (large gray circle) 
and the X-ray point source (tiny gray circle, not to scale) 
at a distance of 1.5 stellar radii, along with surfaces of con- 
stant q (using a velocity law with 7=1 and fo = 0.01). 



The SEI method (Lamers et al. 1987) incorporates an 
exact integration of the radiative transfer equation, yield- 
ing the observed flux F^ at the observed velocity v. 

^=h^S,^^F, (13) 

dT^ 

The optical depth r is parameterised as a function of 

velocity v. 

T f^{ , — for V < vi 

JZ^^^ ~ (14) 

else 




1 



(15) 



(16) 



The parameters ai and 02 are defined as: 

ai=: f — 2 & q;2 = s, 

7 
where t and s parameterize the ionization fraction kj as a 

function of radial distance and velocity: 
Ki(a;) oc a;"^v(a;)"* . (17) 

If m is the dominant stage of ionization and i is the ioniza- 
tion stage of the ion under consideration, then t = m — i. 
If the abundance of the dominant ion in the wind is con- 
stant, one gets 

ai = --2 & ^2=0. (18) 

If wi = 1, the ion is present up to the regions where the 
wind has reached its terminal velocity. 

The intrinsic line profile resembles a GauBian, with 
isotropic broadening a^ due to the thermal and turbu- 



lent motions — hereafter referred to as the "turbulence" 
(Lamers et al. 1987): 

$(Aw) = ^^exp- ( — ) . (19) 

^/TTdy \(Jv J 

This introduces multiple resonance zones in the wind, 
which results in broader and stronger absorption effec- 
tively moving the emission peak towards longer wave- 
lengths (in agreement with P-Cygni line obserservations). 

In the case of doublets (most of the important UV 
resonance lines are doublets), depending on the doublet 
separation the source function of the red component in- 
cludes a contribution from the blue component. Except 
for this coupling, the principle of the SEI code is the same 
as for singlets (Lamers et al. 1987). However, since we 
now include a Stromgren zone we can no longer use a one- 
dimensional radial grid to calculate the source function of 
the red component (as is the case for the SEI method). The 
contribution of the source function of the blue component 
to the source function of the red component depends on 
whether the coupled blue component points are situated 
inside the Stromgren zone. Hence the radial profile of the 
source function depends on the position of the Stromgren 
zone. This is solved by replacing the one-dimensional grid 
when calculating the source function of the red component 
by a two-dimensional axi-symmetric grid. 

The SEI method allows the inclusion of a photospheric 
component in the input line spectrum, but we did not use 
this option in the model profiles shown below. 

The accuracy of the calculations depends on the grid 
/g, which samples the intrinsic line profile. The number 
of points used in the grid is /g,max- In the original code 
a choice of /g,max = 11 yields satisfactory results, but if 
a Stromgren zone is included this value should be some 
two orders of magnitude larger to suppress the numerical 
noise. Stronger turbulence demands larger /g.max- 

3.2. Application of the modified SEI model 

To illustrate how the important model parameters affect 
the shape of the spectral line profile, model spectra are 
shown for a hypothetical resonance doublet of a HMXB 
system. The canonical model consists of a HMXB with 
XX — 1.5, wind velocity law with 7 = 1, wo = 0.01 and 
(Tv — 0.2, doublet separation of 1.5, integrated optical 
depth of the blue component of T = 100 (twice that of 
the red component) , dominant ionization stage in the uni- 
formly ionized undisturbed wind (t, s = 0), no additional 
source of emission (e = 0), no photospheric absorption 
profile, and a Stromgren zone of size q — 3.2. 

Figs. 4 to 8 illustrate the effect of varying several of 
these parameters. The top panels present the UV res- 
onance profiles at two orbital phases: at X-ray eclipse 
{(j) = 0, drawn line) a large fraction of the Stromgren 
zone is behind the OB supergiant so that only a small 
part of the observable stellar wind is affected. When the 
X-ray source is located in the absorbing column in front of 
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Fig. 4. Changing the size of the X-ray ionization zone: SEI model profiles for a HMXB with a wind velocity law 
with 7 = 1, Wo = 0.01 and a^ — 0.2, xx ~ 1-5 (wx — 0.34), doublet separation of 1.5, integrated optical depth 
of the blue component of T = 100 (twice that of the red component), dominant ionization stage in the uniformly 
ionized undisturbed wind, and no additional source of emission. The effect of different sizes q of the Stronigren zone 
is illustrated (closed surfaces for q > 2.941). The top panel shows the profile at X-ray eclipse (drawn line, (j) ~ 0) and 
at ^ = 0.5; the bottom panel presents the difference profile. Note that (much) less absorption is observed at </) = 0.5 
when the X-ray source is in the line of sight. 



the OB supergiant (0 = 0.5), the blue-shifted absorption 
trough is reduced in strength. Note that also the emis- 
sion peak of the P-Cygni profile is affected, though less 
pronounced because the wind volume contributing to the 



P-Cygni emission is much larger in comparison to the vol- 
ume of the absorbing column. The bottom panels show 
the difference spectra {(j) = 0.5 minus </> = 0). 
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Fig. 5. The effect of turbulence: same as Fig. 4, but for q = 3.2 and different levels of wind turbulence ct. 



2 4 

V 




Fig. 6. From weak to strong: same as Fig. 4, but for q = 3.2 and illustrating the effect of increased integrated optical 
depth T. 



The size q of the Stromgren zone strongly affects the larger Stromgren zone enhances the HM-effect especially 

spectral line variability (Fig. 4). Small Stromgren zones near velocities < —vx, also reducing the absorption in 

(large q) leave little trace in the spectral line profile, apart the blue wing of the line profile, and reduces the emission 

from some diminished absorption at = 0.5 at velocities at positive velocities at (f) = (notably in the red com- 

between about —vx and (i.e. not centred at —vx)- A ponent of the doublet). Saturation is not maintained for 
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Fig. 7. Same as Fig. 4, but for q = 3.2 and illustrating the effect of different parameterisations of the optical depth 
(via ai) and the velocity law (via 7). 
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Fig. 8. Same as Fig. 4, but for q = 3.2 and illustrating the effect of different distances xx of the X-ray source 



Stromgren zones that cover a significant fraction of the orbital phases, thereby reducing the contrast in line pro- 
hemisphere in front of the star. However, open Stromgren file shape between transit and eclipse of the X-ray source, 
zones leave only a small area behind the primary unaf- The appearance of the HM-effcet is most prominent for 
fected (in the most extreme case this is a shadow wind), q^^criticai- The definition of q allows a Stromgren zone 
and hence cause diminished absorption and emission at all extending into the X-ray shadow behind the supergiant; 
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obviously, this can only happen if scattering of X-rays by 
ions in the stellar wind is important. 

The presence of turbulence in the wind causes the line 
profile to broaden, which results in less pronounced emis- 
sion (Fig. 5). The HM-effect is also smoothed in velocity 
space. The amplitude of variability at positive velocities 
and between and — Wqo decreases, whilst the amplitude of 
variability in the blue wing {v < — Woo) increases. The blue 
wing variability becomes dominant in the case of strong 
turbulence. In the absence of turbulence it is difficult to 
hide the HM-effect unless the Stromgren zone is much 
smaller than the projected stellar disk. 

Greater optical depth enhances the extent to which the 
absorption is maintained when the Stromgren zone moves 
through the line of sight (Fig. 6). Note that for heavily 
saturated lines with turbulence, as for unsaturated lines, 
it is not straightforward to determine the exact value of 
Woo, and careful modelling is required. 

The ionization stage of the ion, relative to the domi- 
nant ionization stage in the (undisturbed) wind is of im- 
portance for the parameterisation and normalisation of 
the optical depth. If the ion corresponds to an ionization 
stage that is one level below the dominant ionization stage, 
then Qfi — —2 instead of —1 in the canonical example (see 
Eq. 16). For the same integrated optical depth, this would 
yield unsaturated and rather triangular shaped absorp- 
tion profiles with little emission. In the opposite case, if 
the ion corresponds to one level above the dominant ion- 
ization stage, then ai = 0. This yields a heavily saturated 
line profile with stronger emission and somewhat enhanced 
HM-effect in the absorption part between —vx and — Woo, 
compared to the canonical case. From the normalisation 
of the optical depth (Eqs. 14 & 15) it follows that t/T at 
w == -Woo (== vi) is 0.010, 0.22, or 1.0 for ai == -2, -1 and 
0, respectively. Scahng the total optical depth T to 2200 
and 22 for ai = —2 and 0, respectively, the line profiles 
and appearance of the HM-effect (Fig. 7) are similar to 
the canonical case. It may therefore be difficult to distin- 
guish between different values of ai , and prior knowledge 
of the dominant ionization species is required. 

The parameterisation of the optical depth involves the 
velocity law, and hence the parameter 7. Changing 7 from 
1 to 0.5 requires rescaling to q — 1.87, ai = and T = 22 
(Eqs. 4 & 15). The resulting variability is not exactly a 
scaled version of the canonical case (Fig. 7), because q 
depends on both v and x. Most notably, the blue wing HM- 
effect is stronger and saturation more difficult to maintain. 

A different distance xx of the X-ray source implies 
that different parts of the wind are probed. Changing to 
Xx — 1-2 and 2.0 requires rescahng to q ^ 6.2 and 2.15, 
respectively. Increasing the distance of the X-ray source 
enhances the HM-effect in the blue wing at the expense 
of the HM-effect in the absorption part between —vx and 
— Uoo) and also makes it more difficult for a saturated line 
to remain saturated (Fig. 8). If the X-ray source is located 



Table 3. SEI model parameters for fits to the observed 
profiles and variability of UV resonance lines in the spectra 
of HD77581/Vela X-1 and HD153919/4U1700-37. 



Line T ai 


02 


a 


t-oo 


q 


HDll 581 /Vela X-1: 










N v 3 1 


-1 


0.45 


600 


2.9 


Si IV 300 -1 





0.45 


600 


2.9 


C IV 200 





0.45 


600 


2.9 


Al III 20 -1 





0.45 


600 


2.9 


HD 15391 9/4 Ul 700- 37: 










N v 30 





0.15 


1700 


>4 


Si IV 2 X 10* -2 





0.15 


1700 


>4 


C IV 5000 -2 





0.15 


1700 


>4 



close to the primary, then the HM-effect in the emission 
at positive velocities is less pronounced. 

4. Model fits to the ultraviolet resonance lines of 
HD77581/Vela X-1 & HD153919/4U1700-37 

The modified SEI code is used to fit the line profiles at two 
orbital phases (</) = 0, X-ray eclipse, and <}> = 0.5) of N v, 
Si IV, C IV and Al III resonance lines in HD77581/Vela 
X-1 and N V, Si IV and C iv in HD153919/4U1700-37. 
The models are constrained by the following conditions: 
(i) the values for 7, v^a, c and q are the same for all lines 
in a given HMXB; (ii) the value for xx follows from the 
orbital parameters (Table 1); (iii) the dominant ioniza- 
tion stage is estimated from the spectral type (Lamers et 
al. 1999). In HD77581 N2+, 81^+ and C^+ are expected 
to dominate, while in HD153919 probably N'^^, Si'*+ and 
C^+ dominate. AP+ is assumed to dominate in HD77581, 
as the observed Al III line is very strong. The ionization is 
uniform throughout the undisturbed wind, i.e. a2 = (ex- 
cept N''+ in HD77581). Note that the N v hue suffers from 
strong absorption by the wing of Ly-a. A sudden jump in 
intensity near the maximum of the Si iv and C iv emission 
in HD153919/4U1700-37 is due to numerical difficuhies. 
These lines were calculated by sampling the intrinsic line 
profile over ±5cr with /g.max = 3333 (instead of ±3cr with 
-^g,max = 1111), suppressing the spikey feature by half. 
The resulting line profiles are shown in Figs. 9, 10 & 14 
and the model parameters are listed in Table 3. 

The orbital modulation of the strong wind lines in 
HD77581 (with the exception of N v) is mainly due to the 
HM-effect. Some detailed changes are probably related to 
the accretion flow in the system (Sadakane et al. 1985; 
Kaper et al. 1994); since these structures are not included 
in the model, they cannot be reproduced. The modified 
SEI model can naturally explain the (near) absence of 
orbital modulation in the resonance lines of HD153919, 
which is one of the main motivations for this study. 

11. HD77581/Vela X-1 
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HD77581/Vela X-1 




Fig. 9. SEI models that approximately fit the observed line profiles and variability of Si iv, C iv and Al in in the UV 
spectrum of HD77581/Vela X-1. Parameters are summarised in Table 3, and v = 1 corresponds to 600 km s^^. 



4.1.1. Strong wind lines 

Adopting the ionization structure of the stellar wind as ex- 
pected for single stars of the same effective temperature as 
HD77581, and adopting Woo = 600 km s^^, the shape and 
variability of the Si iv, C iv and Al ill (^S-^P° resonance 
doublet: restwavelengths 1854.716 & 1862.790, separation 
1302 km s^^) line profiles are reproduced rather well (Fig. 
9, compare with Fig. 11). 

The N V line profile (Fig. 10, calculated with Ig^max = 
5555), however, cannot be modelled using a normal ioniza- 
tion balance that would yield an extremely profound HM- 
effect, especially at the inner part of the absorption (Fig. 
10a), which is not observed (Fig. 11). Instead, the obser- 
vations seem to favour the degree of ionization to increase 
with distance to the primary: a2 = — 1 gives an acceptable 
fit, suppressing the orbital modulation of the absorption 
trough, while preserving some of the emission (Fig. 10b). 
The line profile and its modulation are very similar for 
«! = or — 1 (Fig. 10c), only requiring a slightly larger 
integrated optical depth (T). If aiJj — 2, however, much of 
the absorption arises from material at extreme velocities, 
causing the absorption profile to extend up to large pos- 
itive velocities beyond the severely diminished emission. 
The absorption observed longward of the N v line can be 
attributed to a C ill line. 

The slope of the blue absorption wing (especially in 
Si IV and C iv, see Fig. 2) indicates that the level of sat- 
uration (deeper and at velocities slightly more negative 



than —Woo) is not reproduced well by the SEI turbulence 
recipe (Eq. 19). The shocked wind structure produced by 
numerical simulations (Owocki 1994), however, seem to 
predict more material at extreme velocities (rather than 
gauBian turbulence). Furthermore, the turbulence may not 
be uniform throughout the wind, as it is assumed here. 



4.1.2. Variations in addition to the HM-effect: evidence 
for a photo-ionization wake 

The main discrepancy between the model and observa- 
tions of HD77581/VelaX-l is the -400 to km s^^ region, 
where the observed intensity is often lower a,t (f) = 0.5 than 
at (/) = 0, contrary to what the HM-effect predicts. This 
may be due to the presence of a photo-ionization wake 
that enhances absorption at phases (f>^0.5 (Kaper et al. 
1994). The photo-ionization wake, as observed in strong 
optical lines like the hydrogen Balmer lines and He i lines 
in the spectrum of HD77581, is situated close to the star 
at the trailing border of the ionization zone. 

This effect is especially pronounced in N v (Fig. 11), 
possibly as a result of a change in ionization degree within 
the photo-ionization wake. Enhanced absorption is seen 
between -|-200 and —600 km s~^ (= — Voo), with its 
maximum evolving towards more negative velocities from 
4> = 0.25 to (j) = 0.75. Absorption at cj) = 0.25 (or enhanced 
emission at </) = 0) extends up to -}-600 km s~^. 
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HD7?581/Vela X-1 




Fig. 10. Same as Fig. 9, but for the N v line in the UV speetrum of HD77581/Vela X-1 (see also Table 3). 



Similar enhanced absorption is seen in Si iv, C iv and 
Al III (Fig. 11) at -200 and -300 to -400 km 8"^ around 
(f) = 0.5 and 0.75, respectively, and perhaps already around 
4> — 0.25. The enhanced absorption due to the photo- 
ionization wake may have reduced the amplitude of the 
HM-effect seen between —600 and —400 km s~^. 

Several weaker UV lines also show orbital modulation. 
The 2po_2D resonance doublet of C ii at 1334.5323 and 
1335.7077 A (Fig. 12) is sufficiently strong to show the 
HM-effect between and about +700 km s~^ and between 
about —400 and —600 km s^^, but slightly enhanced ab- 
sorption at —300 km s~^ around (j> — 0.75 may be at- 
tributed to the photo-ionization wake. 

The strongly variable 3p0_3p niultiplet UV 4 of Si in 
(Fig. 12) consists of 6 components, of which 2 are blended: 
the pattern of variability is repeated five-fold. The absorp- 



tion is diminished at —400 km s at 



0.5 due to the 



HM-effect, but enhanced at —300 km s^^ at = 0.75 
due to the photo-ionization wake. In Fig. 13 a lightcurve 
is plotted for the spectral region between 1292 and 1303 
A — corresponding to —2600 and —70 km s^^, respec- 
tively. Absorption is minimal between (f) = 0.4 and 0.5, 
and strongest between = 0.7 and 0.8. The blending of 
the Si III components makes it difficult to detect variabil- 
ity at —600 km s~^, although there is a hint for dimin- 
ished absorption at = 0.5 at this velocity in the blue- 
most component at 1294.543 A. The ^S-^pO (multiplet 
UV 2) resonance singlet of Si in at 1206.51 A could not 
be studied, due to the severe interstellar extinction and 



the strong Lya geocoronal line at 1215.34 A in the neigh- 
bouring echelle order. The subordinate singlets of ^P°-^S 
multiplets UV 9 and UV 10 of Si in at 1417.24 and 1312.59 
A, respectively, do not show variability, although the lines 
are clearly present in the spectrum — particularly the UV 
9 singlet. The 2p0-2s resonance doublet of Si ii at 1304.372 
A (-1-240 km s~^ in Fig. 12) does not show any variabihty 
either, as it is of interstellar origin. 

The spectral region around 1900 A is dominated by nu- 
merous lines of Fe in, of which the a^S-z^P° multiplet 34 is 
the strongest. It consists of 3 widely separated components 
at 1895.456, 1914.056 and 1926.304 A (Fig. 12). They all 
clearly show the diminished absorption at —400 and —600 
km s~^ around (\> — 0.5 due to the HM-effect, and en- 
hanced absorption at —300 km s~^ around = 0.75 due 
to the photo-ionization wake, exactly as observed in the 
Al in resonance doublet. Some confusion arises from the 
presence of 2 lines of the a^G-z^H° multiplet 51 of Fe ill 
at 1922.789 and 1915.083 A, corresponding to -550 and 
— 1750 km s^^ with respect to 1926.304 A, respectively. 
Multiplet 51 shows the same behaviour of the absorption 
near —300 and —400 km s""'^ as multiplet 34, but not the 
—600 km s^^ absorption. 

In conclusion, the orbital modulation of both the 
strong wind lines and other lines in the UV spectrum of 
HD77581 indicates the presence of two absorption compo- 
nents that are not included in the SEI modelling. The most 
prominent of these can be explained by a photo-ionization 
wake, causing additional absorption in the line-of-sight at 



van Loon et al.: Modelling UV resonance lines in HMXBs 



15 



Vela X-1 KV {XglPA'd.WA] 

■ '■■■■ .. !. r' V 



"^1. 



f: 



\i 



:M 



V 



U,,/' 



r/ V .J^V- 









^^VqAoi 



,^. 



I ' '= - ^ ° ^ s a 



¥'*^f \ -r^'^r 



yoD -:flao 



Vela X-1 KilV (At,t40E.77) 



vn 



'*, h 






I '"^'/^'/fM 



1 ^ / 



^^jkh.4x^"^w^h 












f^fff VS'vV^ 






; ■ — h^ 

l..:^.-^...-..M..Ji f 
I imI 



™>^ D. 






Vela X-1 CW' (,\„1550.774) 



V 



./^a;, 



V \ 



-| ^ ET. ■■ ^0~] 






I ^.s - *° 



Jl^fi.,Sf::\ 



y\ 



'X7^~-' "^VSl 






H 



I i^,^-. -r 



v/VHvvJ.:™ 



■'I^'l' 



f,f^ii^:hiM,„A^ 



aaw -imG 



Vela X-1 Alill (XolS63.79j 



'V\. 



/v.'y\jv\,^ 



A 






\\,A 



rf -''V: 









' ' ' I 

CmZ*£.|- 



_ ' l " 



'!;J- " yVO n:, 



KUW -IDDQ 



, (Km/,] - 



Fig. 11. Phase cj) ^ (solid) and (p — 0.5 (dotted) spectra for the N v, Si iv, C iv and Al ill resonance lines in 
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Fig. 12. Same as Fig. 11, but for C ii, Si iii and Fe iii. 
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HD153919/4U 1700-37 




Fig. 14. SEI models that approximately fit the observed line profiles and variability of N v, Si iv and C iv in the UV 
spectrum of HD153919/4U1700— 37. Parameters are summarised in Table 3, and v — 1 corresponds to 1700 km s~^. 



V ^ —200 km s"'^ around (j) — 0.5 and at w ~ —300 km s^'^ 
around (p = 0.75. In addition, at u ^ —600 km s~^, i.e. the 
terminal velocity of the wind of HD77581, more absorption 
is present outside the Stromgren zone than the SEI model 
reproduces. This might reflect the stellar wind structure 
at larger distances from the star: the velocity may be lower 
and the density higher after the wind has passed through 
the Stromgren zone compared to an undisturbed wind. 



IID153919/4U1700-37. Absorption components that are 
seen both in the strong wind lines and additional weaker 
lines in the UV spectrum of HD77581/Vela X-1 and that 
could not be reproduced by the adapted SEI model can be 
explained by a photo-ionization wake. The orbital modu- 
lation of the resonance lines in the UV spectra of the five 
HMXBs studied here shows a clear trend of the size of the 
Stromgren zone to increase with higher X-ray luminosity. 



4.2. HD153919/4U1700-37 

The shape and (lack of) orbital modulation of the line pro- 
files in HD153919/4U1700-37 are reproduced well (Fig. 
14). The exact value for q is hard to determine because of 
the (near) absence of the HM-effcct due to a combination 
of strong absorption and the presence of turbulence, yet it 
is clear that g < 4 would definitely yield a too strong HM- 
effect. Also, the HM-effect is predicted to appear strongest 
at velocities between — Uoo and 0, implying that any vari- 
ations in the blue absorption wing (at \v\ > Uoo) without 
accompanying HM-effect at smaller velocities must be due 
to some other mechanism, e.g. Raman-scattered far-UV 
emission lines (Kaper et al. 1990). 

5. Discussion 

The adapted SEI model reproduces the HM-effect in 
HD77581/Vela X-1, and naturally explains the lack of 
any (clear) HM-effect observed in the dense wind system 



5.1. Terminal velocities and turbulence 

The terminal velocity Va^ may be estimated from a com- 
parison of an appropriate model line profile with the ob- 
served line profile. It is best derived from moderately sat- 
urated profiles, in which case the maximum absorption in 
the (f> — 0.5 line profile occurs at the terminal velocity. 

The terminal velocity for HD77581/Vela X-1 is thus 
estimated to be Uoo = 600 km s^^, much slower than pre- 
viously reported (1105 km s^^ Prinja et al. 1990). For 
HD153919/4U1700-37 we estimate Woo = 1700 km s^^ 
(Prinja et al. 1990: 1820 km s^^). These values should be 
correct to within 100 km s^^. The low resolution lUE spec- 
tra show that the wind in HDE226868/Cyg X-1 is faster 
(marginally resolved: 1000 < Voo < 1500 km s~^) than in 
Sk-Ph/LMC X-4 and Sk 160/SMC X-1 (unresolved: ^600 
km s^^, see also Hammerschlag-Hensberge et al. 1984 for 
Sk 160/SMC X-1). HST/STIS spectra of the N v. Si iv 
and C IV lines in Sk-Ph/LMC X-4 around (f) — suggest 
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Table 4. Ratios of terminal velocity Woo and escape veloc- 
ity Wesc, and turbulent velocities ay of the primary stars 
in the HMXBs, together with their effective temperatures 
Tcff and F values estimated from Lamers et al. (1995). 



Primary 


Teff (K) 


F 


Vo< 


,/Vcsz 


ay 


HDE226868 


30500 


0.34 


1.8 


—2.8 




Sk-Ph 


36000 


0.16 




^0.8 




Sk 160 


26000 


0.18 




^1.1 




HD77581 


23400 


0.18 




1.3 


0.45 


HD153919 


37200 


0.34 




2.0 


0.15 



Vela X-1 



Silll (X1300) 




Fig. 13. Lightcurve for the integrated flux in the Si in 
resonance lines around 1300 A in HD77581/Vela X-1. The 
dotted line indicates the flux level around d> ~ 0. 



a terminal velocity of ^ 500 km s^^, even though wind 
velocities up to ^ 1200 km s^^ occur (Kaper et al. in 
preparation). 

Single (galactic) 0-type and early-B-type stars are 
known to have a fairly constant ratio of terminal over es- 



cape velocity Woo/fosc ^ 2.5, where 



with 



2GM(1 - F) 



R 



F = 2.658 X 10" 



M 



(20) 



(21) 



for a typical Population I star (e.g. Groenewegen et al. 
1989). This is also true at temperatures between ^ 1.1 
and 2.1 x 10^ K but then t'oo/^'esc ~ 1-3 (Lamers et al. 
1995). The terminal velocities for the primary stars of the 
five HMXBs studied here are, with the possible excep- 
tion of HDE226868, low compared to single stars of the 
same effective temperature (Table 4 & Fig. 15). Especially 
the two magellanic stars have very low terminal velocities, 
which may at least partly be due to their subsolar metal- 
licity and hence lower radiation-force multiplier (Garmany 
& Conti 1985; Prinja 1987; Kudritzki et al. 1987) — which 
is also found for dust-driven winds of Asymptotic Giant 
Branch stars (e.g. van Loon 2000). 



The Stromgren zone could inhibit acceleration of the 
stellar wind material flowing through it, not only re- 
sulting in an ionization wake but possibly also leading 
to a slower terminal velocity of the wind. This scenario 
only works if the Stromgren zone is sufficiently large and 
the orbital period sufficiently short that all wind mate- 
rial leaving the star passes through the Stromgren zone 
before having reached the terminal velocity. The magel- 
lanic systems indeed must have very extended Stromgren 
zones leaving only a small region at the opposite side 
of the primary unaffected by the X-ray source (i.e. a 
shadow wind), and their orbital periods are rather short. 
However, the Stromgren zones in HD77581/Vela X-1 and 
HD153919/4U1700-37 are closed surfaces, and certainly 
do not occupy more than half of the circumstellar space. 
Within half their orbital periods of 4.5 and 1.7 days, re- 
spectively, their stellar winds ought to already have accel- 
erated to velocities approaching the terminal velocity as 
if it were from a single star. 

A viable alternative is that the F factors for single stars 
are not applicable to stars in interacting binaries. HMXBs 
primaries may have already lost a significant fraction of 
their initial mantle mass, causing them to be undermassive 
for their luminosity (Conti 1978; Kaper 2001). This seems 
to be confirmed when comparing the primary star masses 
(Table 1) with Table 3 in Howarth & Prinja (1989). Hence 
we may have under-estimated F and over-estimated Wcgc- 

A turbulence description has been employed to account 
for material at velocities deviating from the bulk flow 
that obeys a velocity law according to standard radiation- 
driven wind theory. The a is the typical deviation in units 
of Woo- To approximately describe both the undisturbed 
wind and X-ray ionized components in the line profiles of 
HD77581/Vela X-1 a very large value for a needs to be 
invoked (cr„ ~ 0.45), whereas for HD153919/4U1700-37 
a much smaller value suffices [a^ ^ 0.15). In absolute 
terms the turbulence in these two HMXBs is very sim- 
ilar, though: Wturb — 270 and 255 km s~^, respectively. 
These values are typical for winds from 0-type stars for 
which Wturb has been found to be largely independent of 
Toff (Groenewegen et al. 1989). As mentioned before, we 
find indications for the deviations from the monotonic ve- 
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Fig. 15. Ratio of terminal over escape velocity for the pri- 
mary stars in the HMXBs studied here (solid dots), com- 
pared to single stars (circles) from Lamers et al. (1995). 
HMXB members have low terminal velocities for their ef- 
fective temperatures. 



locity law to resemble a shocked wind structure rather 
than uniform turbulence. 



5.2. Ionization fractions in the stellar wind 

In principle, the mass-loss rate may be estimated from the 
integrated optical depth T of a resonance line: 

T = hnAuN.v^ (22) 

TOgC 

with /lino the oscillator strength of the transition in terms 
of the harmonic oscillator i^e^ /meC, Aq the rest wavelength 
of the transition, and iVi the column density of the ion: 



M 



ni(r)dr 



(23) 



with n; the number density of the ion. The continuity 
equation yields 



M = ATrr^v{r) (m) 






(24) 
the mean mass of an ion 



with p the mass density, {m) 
in the wind, and Ai the ion abundance by number. The 
mass-loss rate M can now be estimated from observed 
quantities in the following way: 

m) R^v^iP — 



with 



6 /line 



Ao 



(25) 



& 
& 



1-0 = 



Wo > 



(26) 



(1 - 7) for 0<7 < 1 

{vq — l)/lnwo for 7=1 

For our choice of 7 = 1 & wg — 0-01 the numeri- 
cal factor tp — 0.215. Adopting a mean nucleus mass 



(m) = 2.2624 x 10^24 g (Lamers et al. 1999) the prod- 
uct of mass-loss rate M and ion abundance Ai can be 
calculated for the resonance lines in HD77581/Vela X-1 
and HD153919/4U1700-37. 

In practice, however, the limited knowledge of the ion- 
ization balance in the winds of OB supergiants makes the 
derived mass-loss rates highly unreliable. Instead, the ion- 
ization fractions of the ions may be derived when other, 
more reliable estimates for the mass-loss rate are avail- 
able. For resonance lines the excitation fraction of the ion 
that produces the line is unity, and the ion abundance Ai 
is the product of the ionization fraction kj and elemental 
abundance A (by number). Assuming solar abundances 
(Anders & Grevesse 1989), the ionization fraction may be 
derived using Eq. 25. 

Lamers et al. (1999) found kj to depend on both ra- 
diation temperature and density, with the radiation tem- 
perature scaling approximately linearly with Tcs . The em- 
pirical dependencies of Ki on density were contrary to the 
expected ionization balance, however, and it was argued 
that the empirical relations may suffer from selection ef- 
fects. Therefore, we only consider their empirical relations 
between kj and Toff. 

Schroder et al. (in preparation) model the Ha profiles 
of OB supergiants in HMXBs. They derive mass-loss rates 
for HD77581 (M - 1.0 x IQ-e Mq yr^^) and HD153919 
(M ~ l.Ox 10~^ Mq yr^^), which arc comparable to those 
derived for single stars of similar spectral type (Howarth 
& Prinja 1989). Combining their results with the results 
from the SEI models, ionization fractions kj are estimated 
for the N'^+, Si^+ and C^+ (and Al^+) ions in the (undis- 
turbed) stellar winds of HD77581 and HD153919 (Table 
5). These are then compared with the ionization fractions 
estimated from Fig. 3 in Lamers et al. (1999), taking into 
account their lower and upper limits. 

The ion fractions derived from the SEI modelling for 
HD77581/Vela X-1 are in agreement with the predictions 
for single stars by Lamers et al. (1999), except for the N''+ 
abundance which is observed to be an order of magnitude 
higher than predicted. This may be due to either super- 
ionization or nitrogen over- abundance (or both). Auger 
ionization (Cassinelli & Olson 1979) is sometimes invoked 
to explain strong N v resonance lines. Unless the density 
becomes very high. Auger ionization increases with the ve- 
locity in the wind, and may originate in the high-velocity 
extrema of a shocked wind. This might mimic a moderate 
increase of ionization fraction with distance, as required 
to reproduce the observed line profile and variability of 
the N V hue in HD77581/Vela X-1. On the other hand, 
nitrogen over-abundance at the surface of HD77581 may 
have resulted from (i) the transfer in the past of nitrogen- 
enriched material from the progenitor of Vela X-1 onto 
HD77581, or (ii) strong mass loss exposing deeper layers 
mixed with the products of nuclear burning. Kaper et al. 
(1993) note that HD77581 might be a BN star. 
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Table 5. Ionization fractions kj for several ions in the stellar winds of HD77581/Vela X-1 and HD153919/4U1700— 37 
as derived from SEI model parameters, adopting mass-loss rates M from Ha line profile modelling (Schroder et al. 
in preparation) and using solar element abundances A. Ionization fractions from work on single stars (Lamers et al. 
1999) are given for comparison. 



Line Aq (A) fnnc 


T 


A 


Mhc (Mo yr"i) 


Ki 


'^i.Lanicrs 


HD77581/Vela X-1: 

N V 1238.821 0.152 


3 


1.1 X 10"'' 


1 X 10"" 


5 X 10"" 


3 X 10"^ 


Si IV 1393.755 0.528 


300 


3.5 X 10"' 


1 X 10"" 


4 X 10"2 


1 X 10"" 


C IV 1548.20 0.194 


200 


3.6 X lO""* 


1 X 10"" 


6 X 10"^ 


1 X 10"" 


Al III 1854.716 0.560 


20 


3.0 X 10"" 


1 X 10"" 


0.3 


< 1 


HDl 5391 9/4 Ul 700- 37: 
N V 1238.821 0.152 


30 


1.1 X 10"'' 


1 X lO"'' 


2 X 10"-^ 


3 X 10"^ 


Si IV 1393.755 0.528 


2 X 10* 


3.5 X lO""' 


1 X lO"-' 


~ 1 


1 X 10"'^ 


C IV 1548.20 0.194 


5000 


3.6 X 10"" 


1 X lO""' 


6 X 10"2 


3 X 10"-' 



For HD153919/4U1700-37 the opposite is found: the 
observed ionization fraction of N*+ agrees very well with 
the predictions, whereas the ionization fractions of 81"^+ 
and C^^ are observed to be (much) higher than predicted. 
Perhaps these are the dominant ionization states for sili- 
con and nitrogen in the stellar wind of HD153919, rather 
than Si'*+ and C*^ . Still, an ionization fraction k; ^ 1 for 
Si'^"'" is unrealistic for a multi-level atom, and hence the 
SEI model must have over-estimated the integrated opti- 
cal depth of the Si iv line in HD153919/4U1700-37. If 
indeed the degree of ionization in the wind of HD153919 
is lower than that predicted for single stars, this would 
mean that — like in the wind of IID77581 — N^+ is in 
fact overabundant in the wind of IID153919. 

5. 3. Size of the Stromgren zone 

The size of the Stromgren zone, indicated by a particu- 
lar value of the parameter q, is related to the ionization 
parameter ^ as 
qLx 



? 



nxa^ 



(27) 



where a is the distance between the centres of the primary 
and the X-ray source. The ionization parameter can be in- 
terpreted as a measure for the number of X-ray photons 
per particle. The ionization balance in the stellar wind is 
mostly affected by soft X-ray photons. In common with 
previous models of the X-ray ionization of stellar winds 
in HMXBs, we assume the dominant source of opacity 
for soft X-rays is oxygen (Hatchett & McCray 1977). Ma- 
sai (1984) has objected that this assumption may not be 
valid if He Ilis present, but incorporating the effects of a 
He ii/He iii ionization front would be beyond the scope of 
this paper. As a consequence of our assumption, the sharp 
boundaries between the X-ray ionized and the undisturbed 
stellar wind coincide for ions like C*^"*", Si'^"'" and N^+ (cf. 
Kallman & McCray 1982; McCray et al. 1984). Hence the 
edge of the Stromgren zone is given by a particular value 
for ^ that corresponds to the boundary at which oxygen 
is being completely ionized by the X-ray source, and de- 



pends on the shape of the X-ray spectrum and the oxygen 
abundance (Hatchett et al. 1976). 

Accretion of matter onto a star moving through a 
medium was first described by Bondi & Hoyle (1944). 
Their concept was applied to HMXBs by Davidson & Os- 
triker (1973). In an HMXB the compact object has a ve- 
locity Wroi relative to the stellar wind flow: 

^'rcl == f'wind + (bx -V*\ ^Wrot)^ (28) 

a 
with stellar wind velocity fwind near the X-ray source 

(which may be lower than implied by the undisturbed 

wind velocity law; see Section 5.1), orbital velocity vx of 

the X-ray source, and orbital velocity w^, stellar radius R-i, 

and rotation velocity Wrot of the primary. Kaper (1998) 

finds that generally Vj-d ^ ^fwind- Stellar wind matter is 

accreted onto the compact object if it approaches within 

an accretion radius approximately given by 

rel 

with Mx the mass of the compact object. The X-ray lu- 
minosity due to the release of gravitation energy of the 
matter being accreted onto the surface of the compact ob- 
ject with radius i?x is 
GMx 



Lx = TrCr^ccf'roiPx- 



Rx 



(30) 



where px is the mass density of the stellar wind near the 
compact object, and ^ is an efficiency parameter (^ ^ 0.1 
for accretion onto a neutron star). Hence we obtain 

_ ^Rxa^ f Vic\\ .„ ^ 

'^~ 47r({m)\GMx) ^ ' 

The size of the Stromgren zone depends on: (1) the di- 
mensions of the compact object via _Rx, Mx and ^ (by 
the shape of the X-ray spectrum); (2) the orbit via a and 
■yroi (by the orbital velocity); (3) the wind flow via Uroi; 
(4) and the abundances in the wind via (to) and ^ (by the 
oxygen abundance) . The size of the Stromgren zone does 
not, however, explicitly depend on the mass-loss rate nor 
on the X-ray luminosity. 

From the time sequence of the integrated line vari- 
ability (Figs. 1 & 2) and the modelling of the line profile 
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Table 6. Critical, expected and observed sizes of the 
Stromgren zones in the five HMXBs (the smaller q, the 
larger the Stromgren zone). 



HMXB 



ycritical 



Qcxpcctcd yobscrvcd 



HDE226868/Cyg X-1 


1.7 


1.7—5.0 


;$i.7 


Sk-Ph/LMC X-4 


2.7 


<1.6— 4.5 


<2.7 


Sk 160/SMC X-1 


2.9 


<6.3— 17 


<2.9 


HD77581/Vela X-1 


2.7 


10—24 


2.9 


HD153919/4U1700-37 


2.0 


177—215 


>4.0 



variability it is already clear that the Stromgren zone is 
largest for Sk-Ph/LMC X-4 and especially Sk 160/SMC 
X-1: their zones must extend far beyond the primary as 
seen from the X-ray source, and the stellar wind is undis- 
turbed only in a small shadow region behind the primary. 
The Stromgren zone is smaller but still very extended 
for HDE226868/Cyg X-1 where the variability due to the 
HM-effect is continuous over the orbit. HD77581/Vela X- 
1 has an again smaller and now closed Stromgren zone, 
but it still occupies a considerable fraction of the wind 
volume. The (by far) smallest Stromgren zone is found for 
HD153919/4U1700-37. 

The critical and expected sizes of the Stromgren zones 
for the five HMXBs are listed in Table 6. The (/critical is 
calculated assuming a velocity law according to Eq. (2) 
with Vq — 0.01 and 7 = 1. The ^expected is calculated 
assuming ^ = 10'' (Hatchett & McCray 1977) for a typical 
X-ray spectrum, C — 0.1, and Rx — 10 km for Mx = 
1.4 Mq and i?xocMx. The range in values results from 
either assuming co-rotation of the primary with the orbit 
or no rotation of the primary (and the uncertainty about 
Voo for Cyg X-1, LMC X-4 and SMC X-1). The observed 
sizes of the Stromgren zone are in reasonable agreement 
with the expected sizes that, in general, seem to have been 
somewhat under-estimated. This could easily be solved if 
for instance the accretion efficiency C were twice as large. 
In particular, C increases with larger mass and smaller 
radius of the compact object (Shakura & Sunyaev 1973), 
which would yield an accretion efhciency for the black-hole 
candidate Cyg X-1 significantly larger than 10%. 

6. Summary 

An analysis is presented of a large set of lUE 
spectra for five HMXBs: HDE226868/Cyg X-1, Sk- 
Ph/LMC X-4, Sk 160/SMC X-1, HD77581/Vela X-1 
and HD153919/4U1700-37. We compared their spectra 
and variability, and adapted the SEI radiation trans- 
fer code for modelling the variability of resonance lines 
in HMXBs. From model fits for HD77581/Vela X-1 and 
HD153919/4U1700— 37 we derive terminal velocities and 
ionization fractions of their stellar winds, as well as the 
sizes of the Stromgren zone created in the wind by the 
X-ray source. For three other HMXBs (HDE226868/Cyg 
X-1, Sk-Ph/LMC X-4 and Sk 160/SMC X-1) only rough 



estimates are derived for the terminal velocity and size of 
the Stromgren zone. 

The adapted SEI model reproduces the Si iv, C iv 
and Al ill resonance line profiles and their orbital mod- 
ulation in HD77581/Vela X-1, with additional absorption 
components attributed to a photo-ionization wake. The 
N V resonance line in HD77581/Vela X-1 is more diffi- 
cult to model, and it seems to be severely affected by the 
photo-ionization wake. The N v, Si iv and C iv line pro- 
files in HD153919/4U1700-37 are well reproduced by the 
adapted SEI model; their lack of orbital modulation is due 
to the dense stellar wind and low X-ray luminosity. 

OB supergiants in HMXBs have lower terminal veloc- 
ities than single stars of similar spectral type, possibly 
due to a lower effective gravity of the Roche-Lobe fill- 
ing primary. Ionization fractions in the stellar wind of 
HD77581/Vela X-1 agree with predictions for single stars, 
except for an overabundance of N^+ ions in the wind of 
HD77581. This may be due to super- ionization possibly 
resulting from X-ray emission from shocks at the base of 
the wind, or nitrogen overabundance due to mass-transfer 
from the progenitor of Vela X-1 or extensive mass loss 
from HD77581 itself. The N''^ abundance in the stellar 
wind of HD153919/4U1700— 37 agrees very well with pre- 
dictions for single. However, the ionization fractions of the 
other ions (Si^^ and C^^) are much lower than predicted, 
suggesting a generally lower ionization of the stellar wind 
of HD153919 than what is predicted for single stars. This 
would then infer an N'*+ ouer-abundance in the stellar 
wind of HD153919. The sizes of the Stromgren zones are in 
fair agreement with the expectations from standard Bondi 
& Hoyle accretion, and they are larger for more luminous 
X-ray sources. 
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Appendix A: X-ray eclipse spectra 

Averages of spectra near = (eclipse of the X-ray 
source) are a fair representation of the "undisturbed" stel- 
lar wind. The low-resolution spectra of HDE226868/Cyg 
X-1, Sk 160/SMC X-1, and Sk-Ph/LMC X-4 aU contain 
the strong geocoronal Lyman-a emission line at 1215.34 
A, and the resonance lines of C iv near 1550 A and 
Si IV near 1400 A (Fig. 1). These are much stronger in 
HDE226868/Cyg X-1 than in the two Magellanic Cloud 
members that both have a lower metalicity than their 
counterparts in the Milky Way. 

Also immediately clear is the difference between the 
continua of HDE226868/Cyg X-1 and the two Magel- 
lanic Cloud members. Due to the high extinction towards 



22 



van Loon et al.: Modelling UV resonance lines in HMXBs 



Cyg X— 1 = average 




1200 



1400 



1600 



SMC X-1 = average 



1800 
A (A) 




LMC X-4 = average 



1800 
A (A) 




1200 



1400 



1600 



1800 
A (A) 



Fig. A.l. X-ray eclipse spectra of HDE226868/Cyg X-1, 
Sk 160/SMC X-1, and Sk-Ph/LMC X-4. 



HDE226868/Cyg X-1, its continuum is tilted both at the 
1200 A side (dust and hydrogen) and at the 2100 A side 
(grains or large molecules). The Magellanic Cloud mem- 
bers are much less affected by interstellar extinction. 

The high resolution spectra have not been flux cali- 
brated in an absolute sense, therefore the shape of the con- 
tinuum in Figs. A2 and A3 does not represent the shape 
of the continuum as it entered the telescope. 

The geocoronal Ly-a Hne and the resonance lines of 
N V, Si IV, C IV and Al ill are the most prominent features 
in the spectrum of HD77581/Vela X-1 (Fig. A2, top). N v 
and Al iii are very strong for its spectral type. Notice the 
rich Fe ill spectrum at wavelengths longer than ^ 1880 A, 
the Fe iv spectrum between ~ 1480 A and 1880 A, and 
the lack of Fe v and Fe vi at wavelengths ^1480 A. 



In the spectrum of HD153919/4U1700-37 (Fig. A2, 
bottom) the resonance lines of N v, Si iv and C iv and 
the subordinate line of N iv at 1718.551 A are very strong. 
The resonance line of Al ill is mainly interstellar. Notice 
that the Fe ill and Fe IV spectra are absent, whilst at 
wavelengths ;$1480 A the spectrum is crowded with lines 
of Fe V and Fe VI. This contrast with HD77581/Vela X- 
1 is due to the higher photospheric temperature of the 
06.5 Iaf+ star HD153919 compared to the BO. 5 lab star 
HD77581. The stronger P-Cygni hues of Si iv and C iv in 
HD153919/4U1700-37 reflect the denser and faster stellar 
wind of HD153919 compared to HD77581. Note that the 
N V emission is stronger in HD77581/Vela X-1. 



Appendix B: Normalisation and continuum 
variability 

Comparison of column densities towards the primary as 
derived from line absorption requires the continuum to 
be normalised to eliminate continuum variability of the 
primary — which is mainly due to tidal deformation of 
the primary, with a possible contribution of X-ray heating 
of the primary by the X-ray source. It must be realised 
that normalising the continuum is in principle incorrect 
when studying line emission: at a certain orbital phase, the 
emission is due to scattered light originating from parts of 
the primary that are not necessarily the part from which 
the observed continuum originates. 

For each object, three wavelength intervals were cho- 
sen that show no intrinsic variability other than a possi- 
ble continuum variability. These are listed in Table Bl. 
To normalise the flux scales of the spectra relative to ea- 
chother, each spectrum k was divided by a normalisation 
factor /fc. This factor was constructed from the spectral 
points {xki}^^i in the three intervals {{Aj,Bj)}^^^, with 
interval j having rij spectral points: 

3 



fk 



E 



i=Ai 



•^ki 



X 



^N 



-1 2^i=Ai 



(B.l) 



In this way the spectra are corrected for continuum vari- 
ability as long as the continuum slope remains constant. 
The normalisation constants may in principle be used 
to construct UV continuum lightcurves, after correcting 
for the change in sensitivity of the detector according to 
the correction table of Bohlin & Grillmair (1988). The cor- 
rection procedure does not discriminate between the two 
resolution modes (Cassatella et al. 1994). Not all spectra 
can be used to construct a lightcurve. For the HMXBs 
in low resolution all spectra available are used except for 
SWP3989 (HDE226868/CygX-l) and SWP1458 and 1459 
(Sk-Ph/LMC X-4) that had a large (factor of 3 to 5) offset 
in flux level. For the HMXBs in high resolution all spectra 
of HD77581/Vela X-1 are used, but SWP1476, 1714 and 
1972 through 5180 (HD153919/4U1700-37) are omitted 
because they were taken through the small aperture caus- 
ing the loss of an unknown fraction of the total light. 
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Fig. B.l. Lightcurves of the UV continuum flux and the UV colours as defined in the text. 
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Table B.l. Intervals {A,B), in A, used for constructing 
the continuum lightcurves. 



HMXB 



Band 1 



Band 2 



Band 3 



Cyg X-1 
LMC X-4 
SMC X-1 
Vela X-1 
4U1700-37 



(1265,1295) 
(1650,1670) 
(1605,1635) 
(1272,1278) 
(1465,1470) 



(1759,1765) 
(1690,1710) 
(1746,1762) 
(1492,1498) 
(1515,1518) 



(1820,1895) 
(1820,1895) 
(1820,1895) 
(1696,1702) 

(1772,1777) 



To study variability in the spectrum tilt, wc define UV 
colours in the following way: 

C/far = Ll-i2 & Un,id^L2'L3 (B.2) 

where Li, L2 and Lg, are the fluxes (normalised to unity) 
for the wavelength bands 1, 2 and 3 in order of increas- 
ing wavelength. Because of the difficulty to find suitable 
wavelength bands of the best quality for each HMXB, the 
bands — and therefore the colours — are defined for each 
HMXB individually. 

The UV continuum lightcurves are displayed in Fig. 
Bl. The iavo hghtcurve is an average of all points within 
the 3 wavelength bands. It is not straightforward to deter- 
mine the errors on the individual points in the lightcurve. 
The photometric accuracy of lUE is reported to be ~6% 
at a 95% confidence level (Bohlin ct al. 1980). This is 
probably a conservative estimate, judging from the small 
degree of scatter. For the error on the colour measure- 
ment we adopt the standard deviation of the distribution 
function of the measured C/mid points, which is a very con- 
servative error estimate, and the same value is assigned to 
the error of the C/far colour. 
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Fig. B.2. Continuum hghtcurves of Sk-Ph/LMC X-4 for 
orbital phase bins as indicated in each of the nine frames. 
The cycle parameters are those of the processing disk: 
00 = JD 2 443 392.52(2); P ^ 30"? 38(3). The dotted co- 
sine waves are fits- by-eye to the data (see text). 



B.l. HDE226868/Cyg X-1 

The hghtcurve of HDE226868/Cyg X-1 (Fig. Bl) has a 
minimum to maximum amplitude of ~17% (the deviating 
point at = 0.281 probably does not contain all the flux 
of the source). This is consistent with the scatter of similar 
amplitude in the UV photometry around 1500 and 1800 
A obtained by Wu et al. (1982) with the Astronomical 
Netherlands Satellite (ANS), but considerably larger than 
the 8% found by Treves et al. (1980) who analysed only 7 
spectra by integrating the entire spectrum between 1250 
and 1900 A. The absence of X-ray eclipses limits the in- 
clination to « ^ 60° (Bolton 1975). Hence the deprojected 
amplitude is even larger, implying a severe deformation of 
HDE226868. 

The Ufar colour shows clear orbital modulation similar 
to the average flux: when the HMXB becomes fainter in 
the UV, the spectrum becomes redder. This may reflect 
a lower photospheric temperature of HDE226868 at the 
side that faces to and away from Cyg X-1, probably a re- 
sult of the tidal deformation of the primary (Hutchings 
1974). HDE226868 is possibly flUing its tidal lobe dur- 
ing peri-astron of the slightly eccentric orbit (e ^ 0.05: 
Bolton 1975). No heating of the photosphere of the pri- 



mary by the X-ray companion is observed, in agreement 
with Stromgren photometry by Hilditch & Hill (1974). 

B.2. Sk-Ph/LMC X-4 

The hghtcurve of Sk-Ph/LMC X-4 (Fig. Bl) has a mini- 
mum to maximum amplitude of ^19%, and is an improve- 
ment over the lightcurve derived from only 15 spectra by 
van dcr Khs ct al. (1982) (see also Vrtilek et al. 1997). 
The variable depth of the 4> = 0.5 minimum was explained 
by a precessing accretion disk (van dcr Klis et al. 1982). 
The precession causes variations in the amount of sur- 
face of Sk-Ph that is exposed to X-ray heating, as well 
as variations in the fraction of Sk-Ph that is eclipsed by 
the disk. Heemskerk & van Paradijs (1989) confirmed the 
presence of an accretion disk that processes with a period 
of 30.36±0.02 days, whilst X-ray observations had already 
revealed a 30.42±0.03 days period (PakuU et al. 1985); we 
adopt P = 30.38 ±0.03 days. Phase 0prec = corresponds 
to the phase when the accretion disk is seen edge-on. 

Lightcurves are created according to the cycle parame- 
ters of the precession of the accretion disk, for bins around 
orbital phases 0orb = 0, 0.25 & 0.75, and 0.5 (Fig. B2). For 
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orbital phases (/)orb ~ no clear variability with the disk 
precession period is seen. For (/)orb ~ 0.5 the luminosity is 
lowest when the accretion disk eclipses the largest part of 
Sk-Ph (0proc — 0.5), and highest when the accretion disk 
is seen edge-on {(j>prcc ~ 0). This is especially clear when 
narrower orbital phase bins arc taken. At (/)orb ^ 0.25 & 
0.75 the luminosity varies in anti-phase with the varia- 
tion at (j)orh = 0.5. Our results agree well with those of 
Heemskerk & van Paradijs (1989) from optical data. The 
amplitude of variability of the UV luminosity at 0orb = 0.5 
is 10% — the same as in the optical. 

The C/far colour (Fig. Bl) is modulated with the orbital 
period. The colour behaviour around orbital phase (jjorb ^ 
0.5 may be dominated by X-ray heating, causing the UV 
spectrum to become hotter close to (/)orb ~ 0.5. The scatter 
in the Ufar colour very close to 0orb = 0.5 may be due to 
variable obscuration by the precessing disk. 



B.3. Sk 160 /SMC X-1 

The lightcurve of Sk 160/SMC X-1 (Fig. Bl) has a 
minimum to maximum amplitude of ^17%. In optical 
lightcurves the (j) = 0.5 minimum is observed to be weaker 
than the = minimum, explained by X-ray heating by 
the bright X-ray source (Hutchings 1974). Van Paradijs 
& Zuiderwijk (1977) and Howarth (1982), however, found 
that X-ray heating is not sufficient and that an emitting 
accretion disk is required, in agreement with the sugges- 
tion that the mass transfer is dominated by Roche-lobe 
overflow (Hutchings et al. 1977). Despite using more spec- 
tra than van dcr Klis et al. (1982), the orbital phase cov- 
erage of our UV lightcurve around = 0.5 is too poor to 
answer the question of X-ray heating and accretion disk. 
An ^ 60 day periodicity in the X-ray characteristics of 
SMC X-1 has been attributed to a precessing accretion 
disk (e.g. Wojdowski et al. 1998). 



B.4. HD77581/Vela X-1 

The lightcurve of HD77581/Vela X-1 (Fig. Bl) has a min- 
imum to maximum amplitude of ^20% and clear orbital 
modulation, although Dupree et al. (1980) could not see 
continuum variability in the 5 low resolution spectra they 
used. The (f> = 0.5 minimum is significantly deeper than 
the (f> = minimum, much alike the visual lightcurve 
(Zuiderwijk et al. 1977) that has a somewhat smaller am- 
plitude of ^15%. Zuiderwijk et al. argue that HD77581 is 
nearly filling its Roche lobe. 

The C/far colour is modulated with the orbital pe- 
riod, becoming redder around (f> = 0.5. This may be ex- 
plained by temperature gradients over the photosphere of 
HD77581. The scatter of the f/far colour near (p = 0.5 
may result from variable X-ray heating of the part of the 
photosphere of HD77581 that faces Vela X-1, due to the 
strongly variable X-ray flux. 



B.5. HDl 5391 9/4 Ul 700- 37 

The hghtcurve of HD153919/4U1700-37 (Fig. Bl) is flat 
within ^4 to 5%, possibly with a marginable maximum at 
(j) = 0.5. This is less than the UV amplitude around 1500 
and 1800 A of 8% found by Hammerschlag-Hensberge & 
Wu (1977) using the ANS. Optical amplitudes are ~4 
to 8% (Hammerschlag-Hensberge & Zuiderwijk 1977; van 
Paradijs et al. 1978). The small amplitude may be due 
to a large mass ratio of the system and therefore small 
distortion of HD153919. The optical lightcurve cannot 
be described by tidal distortion of HDl 539 19 alone (van 
Paradijs et al. 1978). In particular the deepest minimum 
does not occur at = 0.50 but at </) = 0.56. Possible ex- 
planations include extra absorption originating from the 
mass flow in the system. The UV lightcurves presented 
here also suggest a minimum near (j) — 0.6. Hints for the 
signature of a photo-ionization wake in Stromgren pho- 
tometry are presented in Appendix C. The t/far colour 
is bluest at orbital phase = 0.5, which suggests X-ray 
heating of the part of the photosphere of HD153919 facing 
4U1700-37. 



Appendix C: Stromgren uvby/3 photometry of 
HD153919/4U1700-37 

In August 1993 we obtained Stromgren uvby plus broad 
and narrow band H/3 photometry using the 50 cm 
SAT at ESO/La Silla. The 09.5 lab star HD154368 
was measured for comparison. The photometry is pre- 
sented in Tables CI & C2. The magnitude differences 
between HD153919/4U1700-37 and HD154368 are plot- 
ted in Fig. CI, together with the colours mi and ci 
for HD153919/4U1700-37 as described in Sterken et al. 
(1995). The comparison star did not vary with time. The 
standard deviation for each individual measurement is 
typically 0.005 mag in u and y, and 0.004 mag in v and 
b. From the scatter on small orbital phase scales in the 
resulting lightcurves a formal, conservative estimate of 1- 
(7 ^ 0.01 mag is derived. Because of the advantages of 
the multi-channel photometer, the errors in mi and ci are 
smaller than this. The standard deviation of the distribu- 
tion of measured values gives an upper limit to the error 
on each individual measurement of l-a = 0.0065 mag for 
both colours. The /3 parameter is derived from the H/3 
narrow- and wide-band magnitudes. It increases with in- 
creasing absorption strength of the H/3 line. The data give 
an upper limit to the error in /? of 1-cr = 0.0019 mag. 

Our uvby lightcurves of HD153919/4U1700-37 
agree with the lightcurves obtained by Hammerschlag- 
Hensberge & Zuiderwijk (1977). They noticed an orbital 
modulation of the colour toi, with a minimum at ^ ~ 0.2, 
and suggested that this might be caused by variability 
in the emission lines of He ii 4686, C ill 4650 and N III 
4634-4641 that are included in the Stromgren b-filter. Our 
observations indeed confirm a minimum in toi between 
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Fig. C.l. Stromgren uvby photometry, the derived colours 
TOi and ci, and the /3 parameter (see text) for 
HD153919/4U1700-37. The uvbyP are presented differ- 
entially with respect to the comparison star HD 154368. 
Mean levels are represented by dotted lines. 

(j) ~ 0.1 and 0.2. Krzcmihski (1976) did not find any or- 
bital modulation of the He ii 4686 line. Kaper et al. (1994) 
found variability in the He ii 4686 line profile, but this 
alone cannot cause observable variability in the Stromgren 
b-filter. The other two mentioned emission lines may still 
be variable. 

There is a weak indication for a gradual decrease in the 
H/3 absorption (or an increase in emission) towards (j) — 1, 
possibly with a secondary weak minimum around (f> — 0.3 



Table C.l. Heliocentric Reduced Julian Date 
RJDq =JD-2 449 200, u, v, b and y magnitudes 
(with their standard deviation of the last digit (s) within 
parentheses) and P parameter of HD153919/4U1700-37. 
Orbital phases may be calculated from: (po = JD 
2 446161.3400(30); P = 3':'411652(26) (Haberl et al. 
1989). The last row gives the mean values and their 1-cr 
errors. 
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6.743(4) 


6.556(4) 


1.335 


21.5959 


7.155(3) 


7.051(3) 


6.755(3) 


6.566(4) 


1.183 


22.5641 


7.129(3) 


7.025(2) 


6.731(3) 


6.541(3) 


1.096 


23.5721 


7.187(7) 


7.084(5) 


6.784(4) 


6.595(4) 


1.124 


24.5867 


7.132(4) 


7.025(4) 


6.728(2) 


6.539(3) 


1.180 


25.5770 


7.183(7) 


7.089(8) 


6.787(8) 


6.597(12) 


1.158 


26.5891 


7.161(8) 


7.052(3) 


6.754(2) 


6.566(5) 


1.212 


28.5688 


7.141(3) 


7.044(6) 


6.748(7) 


6.560(7) 


1.159 


32.5778 


7.143(6) 


7.042(3) 


6.745(3) 


6.554(3) 


1.236 


Aug'93 


7.146(7) 


7.043(7) 


6.747(7) 


6.559(7) 


1.233(27) 



(Fig. CI). This is exactly what Kaper et al. (1994) found 
in their study of the variability in the H/3 P-Cygni fine 
profile. They attributed this to the presence of a photo- 
ionization wake. 



Appendix D: Covariances method for error and 
variability analysis 

The degree of variability in a particular spectral point over 
a series of spectra can be expressed by the variance in that 
spectral point. The variances are often used to derive error 
estimates on the values of spectral points, and to detect 
intrinsic spectral variability. Usually two adjacent spec- 
tral points do not behave independently from eachother. 
Ignoring the covariability of adjacent spectral points can 
lead to serious under-estimation of the error on a quantity 
derived by integration along the spectrum. A correct error 
analysis therefore involves the calculation of covariances. 
Better than variances, covariances are powerful in detect- 
ing variability on an intermediate spectral scale — i.e. ex- 
ceeding the instrumental profile but considerably smaller 
than the total spectral range. Covariances have also been 
used successfully in proving the reality of spectral features 
with low signal-to-noise (van Loon et al. 1996). 

For a set of N spectra, each consisting of a number of 
spectral values x, the covariance of points a & 6 is: 

2 ^EfcLl XkaXkb - (EfcLl ^ka){Y.k=l ^kb) 
'^ab = 



N{N - 1) 



(D.l) 
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Table C.2. Same as Table CI, but now for the comparison 
star HD154368. 



Table D.l. Intervals {A,B) 
tion. 



in A, used for error estima- 



RJD0 


u (an) 


V (cr„) 


b{at) 


V i^y) 


/3 


09.6355 


7.224(4) 


6.962(2) 


6.600(2) 


6.151(2) 


1.197 


09.6475 


7.224(4) 


6.959(3) 


6.598(2) 


6.150(4) 


1.250 


13.6283 


7.210(3) 


6.946(2) 


6.585(3) 


6.139(2) 


1.214 


13.6405 


7.207(5) 


6.946(3) 


6.584(2) 


6.140(3) 


1.271 


14.6676 


7.188(4) 


6.928(3) 


6.566(3) 


6.124(3) 


1.463 


14.6795 


7.186(7) 


6.929(3) 


6.567(3) 


6.124(3) 


1.566 


15.6101 


7.215(9) 


6.949(5) 


6.587(5) 


6.139(3) 


1.165 


15.6237 


7.207(4) 


6.947(3) 


6.583(2) 


6.138(3) 


1.219 


16.6434 


7.189(11) 


6.927(5) 


6.562(5) 


6.116(5) 


1.334 


16.6560 


7.182(4) 


6.924(4) 


6.560(5) 


6.116(5) 


1.418 


17.6364 


7.183(5) 


6.918(4) 


6.558(4) 


6.112(4) 


1.310 


17.6484 


7.177(6) 


6.920(4) 


6.558(4) 


6.117(4) 


1.384 


18.5696 


7.203(5) 


6.931(4) 


6.575(4) 


6.142(5) 


1.075 


18.5832 


7.209(6) 


6.940(4) 


6.582(4) 


6.148(4) 


1.107 


19.6155 


7.193(33) 


6.933(8) 


6.570(8) 


6.127(6) 


1.232 


19.6283 


7.185(6) 


6.926(6) 


6.567(5) 


6.123(6) 


1.296 


20.6244 


7.199(6) 


6.933(4) 


6.572(4) 


6.126(4) 


1.291 


20.6374 


7.192(6) 


6.933(5) 


6.572(4) 


6.127(5) 


1.368 


21.5897 


7.198(3) 


6.939(2) 


6.575(3) 


6.127(3) 


1.153 


21.6021 


7.206(6) 


6.939(6) 


6.575(5) 


6.127(5) 


1.199 


22.5579 


7.198(4) 


6.935(3) 


6.572(2) 


6.124(2) 


1.074 


22.5707 


7.195(3) 


6.938(3) 


6.574(3) 


6.129(3) 


1.104 


23.5659 


7.204(6) 


6.938(5) 


6.577(4) 


6.129(4) 


1.099 


23.5784 


7.197(4) 


6.939(3) 


6.576(3) 


6.132(4) 


1.134 


24.5792 


7.215(3) 


6.956(2) 


6.594(2) 


6.148(3) 


1.146 


24.5931 


7.213(4) 


6.952(3) 


6.592(2) 


6.145(5) 


1.197 


25.5690 


7.234(9) 


6.972(7) 


6.607(7) 


6.158(7) 


1.124 


25.5837 


7.236(12) 


6.987(15) 


6.622(15) 


6.182(18] 


1.172 


26.5828 


7.241(8) 


6.975(5) 


6.611(4) 


6.161(4) 


1.179 


26.5956 


7.230(6) 


6.970(3) 


6.607(3) 


6.162(5) 


1.231 


28.5619 


7.206(11) 


6.948(4) 


6.584(4) 


6.137(4) 


1.128 


28.5756 


7.204(4) 


6.945(2) 


6.582(2) 


6.137(3) 


1.173 


32.5704 


7.224(6) 


6.969(4) 


6.607(4) 


6.164(6) 


1.196 


32.5859 


7.210(8) 


6.951(6) 


6.589(7) 


6.146(6) 


1.266 


Aug'93 


7.205(3) 


6.944(3) 


6.582(3) 


6.137(3) 


1.228(20) 



Cyg X-1 LMC X-4 SMC X-1 Vela X-1 4U1700-37 



It measures how much and how coherently the values of 
the two points vary from one spectrum to another. For a ~ 
b the covariance reduces to the variance. Integrating the 
spectral values x over a spectral range {A, B) in spectrum 
k yields: 



Xkidi 



(D.2) 



with an error estimate cr/^ given by: 

^l =11 4djd^ (D-3) 

J A J A 

Only if all points within {A, B) behave statistically inde- 
pendently this reduces to the commonly used variance- 
deduced error estimate, because then: 
<7ij = o-iiSij (D.4) 

with Sij the Kronecker delta function. In the opposite ex- 
treme when all points within (A, B) behave in phase: 
(^ij = o-ii (D.5) 



1250,1370 
1420,1510 
1565,1926 



1170,1200 
1280,1380 
1430,1530 
1565,1926 



1250,1380 
1415,1530 
1560,1926 



1153,1200 
1253,1280 
1315,1380 
1415,1530 
1565,1840 



1183,1225 
1420,1517 
1730,1947 



and consequently no increase in signal-to-noise can be ob- 
tained by integrating along the spectrum. The integrated 
covariances J^ erf dj are calculated in spectral regions 
without intrinsic spectral variability (Table Dl) and tab- 
ulated in flux-wavelength space. Errors can then be as- 
signed according to flux and wavelength. This assumes 
that the integrated covariances are a smooth function of 
flux and wavelength, which was proven to be true for lUE 
spectra (Howarth & Smith 1995). 

The covariances cr? decrease with increasing distance 
\j — i\ to the point i. If the integral of the covariances 
for the point i converges sufficiently fast, then the inte- 
gration interval {A, B) in Eq. (D3) may be replaced by a 
smaller interval such that the integral of the covariances 
just reaches convergence. This interval is determined from 
a covariance profile of cr? versus j — z, which represents the 
spectral shape of the instrumental profile. As an example 
the covariance profile and the cumulative covariance are 
shown for Cyg X-1 and Vela X-1 (Fig. Dl). The dashed, 
dotted and dash-dotted lines represent the mean curves 
within the three regions in which the continuum calibra- 
tion factors were determined (from shortest to longest 
wavelength, respectively), whilst the solid line represents 
the entire spectral range — i.e. including wildly variable 
spectral features. Both the covariance and the cumulative 
covariance at zero distance from the point i reduce to the 
variance — i.e. the covariance of point i with itself. At 
larger distances the covariance diminishes to a small os- 
cillation around zero, whilst the initially rapidly growing 
cumulative covariance approximates a constant level. 

For Cyg X-1 the solid line is too high to be captured 
within the frame of the plots, indicating that some parts 
of the spectrum are strongly variable. The width of the 
covariance profile corresponds to the specified spectral re- 
solving power of about 250. The cumulative covariance 
indicates that a reasonable choice for the integration in- 
terval is ±6 A. The sharp peak of the covariance profile 
of Vela X-1 indicates a spectral resolution in accordance 
with the specified spectral resolving power of about 10^. 
Its broad wings, however, prevent the cumulative covari- 
ance from converging within ^16 spectral points (= 1.6 
A) — several times the specified spectral resolution. 

Spectra of the variances and integrated covariances are 
shown for Cyg X-1 and Vela X-1 (Fig. D2) as an exam- 
ple of the detection of intrinsic spectral variability. In the 
low-resolution spectra of Cyg X-1 the spectral features 
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Fig. D.l. Covariance (top) and cumulative covariance (bottom) for Cyg X-1 (left) and Vela X-1 (right) as a function 
of distance along the spectrum with respect to the spectral point i. The drawn, dashed, dotted and dash-dotted hues 
represent the averages for the entire spectrum and the three spectral regions in which the continuum calibration factors 
were determined, respectively, with subsequently larger mean wavelength. The vertical long-dashed line indicates the 
choice of the integration boundary. 



(notably the Si iv and C iv resonance lines near 1400 and 
1550 A, respectively) and their variability are unresolved. 
The calculation of covariances does not improve the de- 
tection sensitivity much. In the high-resolution spectra of 
Vela X-1, however, the integrated covariance spectrum is 
much more sensitive than the variance spectrum in detect- 
ing intrinsic variability. Events that are limited to only one 
or very few adjacent pixels (like cosmics or dead pixels) 
cause a "forest" of sharp peaks in the variance spectrum, 
nearly completely obscuring the resolved but still narrow 
peaks of the intrinsic variability — the "trees" of inter- 
est. By calculating the integrated covariance spectrum the 
forest is suppressed, and at the same time the features in- 
dicating intrinsic variability are enhanced. The variable 
spectral features can be identified easily. 
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